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Abstract 

We propose that the correlation functions of the inhomogeneous eight- 
vertex model in the anti-ferroelectric regime satisfy a system of difference 
equations with respect to the spectral parameters. Solving the simplest 
difference equation we obtain the expression for the spontaneous staggered 
polarization conjectured by Baxter and Kelland. We also discuss a related 
construction of vertex operators on the lattice. 



Introduction 

In ||l|, ^ it was recognized that the correlation functions of the inhomogeneous 
six-vertex model in the anti-ferroelectric regime can be expressed as a trace of 
products of the g-deformed vertex operators. An explicit integral formula is 
given in ||^ . These correlators satisfy a system of g-difference equations , that 
were introduced by Smirnov in the study of form factors in massive integrable 
QFT Q and correspond to the level-0 case of the q-KZ equation of Frenkel and 
Reshetikhin g. (To be precise the equations for the correlators are 'dual' to 
Smirnov's ones, but we do not go into such details here). As remarked in if 
one replaces the trigonometric R matrices appearing here by the elliptic ones, 
the resulting equations are still 'completely integrable'. In this paper we propose 
that the latter are precisely those satisfied by the correlators of the eight- vertex 
model in the anti-ferroelectric regime. 

Our method is based simply on the Yang-Baxter equation and the crossing 
symmetry and, as such, is applicable to more general models. This construction 
also allows one to interpret the g-vertex operator employed in Q ||] as an oper- 
ator that inserts a dislocation (an extra half infinite line) on the lattice. In this 
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formulation the vertex operators generalize straightforwardly to the elliptic case. 
It remains an interesting open problem to give a mathematical construction of 
such operators, along with the relevant elliptic deformation of Lie algebras. 

The text is organized as follows. In Sect. 1, we introduce our notation, and 
formulate the properties of the general correlators including Smirnov's difference 
equations. We give a heuristic argument for derivation of the main statements. 
In Sect. 2 we solve the simplest difference equation and derive the spontaneous 
staggered polarization. We obtain the expression conjectured by Baxter and 
Kelland [^. Sect. 3 is devoted to the construction of vertex operators on the 
lattice. 



1 Correlators and Smirnov's equations 



1.1 The eight-vertex model Consider an infinite square lattice consisting of 
oriented lines, each carrying a spectral parameter varying from line to line. The 
orientation of each line will be shown by an arrow on it. A vertex is a crossing 
of two lines with spectral parameters, say Ci and C2, together with the adjacent 
4 edges belonging to the crossing lines, as shown in (|^). The edges are assigned 
state variables (which we call 'spins') : ei,e2,£i,£2. In the eight- vertex model, 
each spin can take one of the two different values ±. A spin configuration around 
the vertex is an assignment of ± on the four edges. Notice that our notation 
is somewhat different from the standard. In particular, arrows on the lines 
denote the orientations that we assign to them, rather than the spins. There 
are 16 possible vertex configurations. We assign each configuration a Boltzmann 
weight. The set of all Boltzmann weights form the elements of the i?-matrix: 



£1 



Ci 



-S2 = R{Cl/C2)e[, 



(1) 



The matrix i?(C) acts on (g) via the natural basis of 



as 



per site k, it is given explicitly by 



Normalized by the partition fmiction 



RiO = 



<0 



/a{C) d{0\ 

b{c) c(c) 
c(c) &(C) 



— \d{0 







a(C)/ 



where the unnormalized Boltzmann weights a, b, c, d (see (10.4.23-24) and (10.7.9) 
in 0) are given by 



«(C) 



-ipQ{i\)H{i\-iu)<d{iu), 6(C) 
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c(C) = -ipH{iX)e{i\-iu)e{iu), d{C) ^ipH{i\)H{i\-iu)H{iu), 



where we set q = exp(—7rl'/l), x = ea::p(— 7rA/2/), ( = exp{'Ku/2I). We shall 
restrict our discussion to the principal regime 

< g < a; < < 1, 

in which k(() is given by (see (10.8.44), ibid.) 

i^iO = pjx~^K{C^){x'^C'^;q)oo{qx^^C^;q)oo, i ^ q^'^{q]q)lo{q^',q^)oo, 

{x'^z; q, x'^)oc{x'^z~^;q, x*)oo{qz; q, x^)oc{qx'^ z'^; q, x'^)ao 



k{z) 



{x^z ^■,q,x'^)ocix'^z;q,x^)ooiqz '^]q,x^)ao{qx'^z]q,x^)c 



The symbol [z; qi, . . . , qk)oc means n^....,nfc=o(l ~ ' • • • qT)- Notice that in 
the definition of the i?-matrix, k was used to normalize the Boltzmann weights. 
The _R-matrix satisfies 



R{OPR{C')P = l. 

{PR[C^)P)'' = (a" ® l)i?(C/a;)(a" ® 



(2) 



Here Pv ® v' 
component, and 



v' (E) V, (•)*^ means the transpose with respect to the first 

r 

1 



We note that 



i?(x-i) = 
and that Eq.(||) is equivalent to 



£l 








1 — < — 


e 


/ 

1 



-£2 









1 
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Vo 






£l 




— < — 




C2 




£'i 



0^ 



(3) 



Namely, reversing the orientation of a vertical line in ([Tj) gives rise to reversing 
the spins on that line and shifting the spectral parameter Ci ~^ Ci = Ci/x- 
1.2 Correlators and dislocations on the lattice Let us consider an infinite 
square lattice Q with all the vertical lines oriented downward, all the horizontal 
lines to the left. We associate the spectral parameters Q to the vertical lines, 
and $,k to the horizontal lines. As argued in the calculation of correlators 
of arbitrary spins reduces to calculating correlators of the vertical-edge spins 
located on the same row, where by 'row' we mean the set of all vertical edges 
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between two neighboring horizontal hncs. We shall restrict our attention to this 
case. 

In the principal regime we are considering, the Boltzmann weight c domi- 
nates the others. In the low temperature limit: q, x, S,k/Cj ~^ 0; only type c 
vertex-configurations are non-vanishing. Therefore, the spin variables take the 
same value in the NE-SW direction, and alternate in the NW-SE direction. In 
this limit, two spin configurations are possible: the ground state configurations 
(GSC's). At finite but low temperatures, we choose and fix a GSC, and con- 
sider the statistical sum over all configurations which differ by a finite number 
of spins from that GSC. We let « G Z2 label the choice of GSC by the fact that 
in the low temperature limit the spin on some reference edge is frozen to the 
value (—1)*+^. The choice of the reference edge will be given later. Now choose 
a particular row, and consider n successive vertical edges on it. Let C,i, . . . Xn 
be the spectral parameters attached to the corresponding lines, numbered from 
left to right. As the reference edge we take the next left to the one with the 
spectral parameter (^1 . Consider the probability of the spins taking the values 
£1, . . . , e„ e {±}. Baxter's conclusion |^ was that this correlator is independent 
of all other spectral parameters. 

In this work, we consider sums over all spin configurations on various lattices, 
with certain spins kept fixed. On infinite lattices, such configuration sums are 
not well-defined. Only the ratios of such sums, on the same lattice but with 
the different choices of fixed spins, are well-defined in low temperature series 
expansions. When considering probabilities, the normalization is so chosen that 
the sum of all correlators is equal to 1. However, except for the simplest situation 
as (|4|), such a choice does not give solutions to our system of equations. In fact, 
we do not know how to normalize our correlators. We postulate that there exists 
a proper normalization such that they satisfy all the required properties given 
below. 

In order to obtain expressions that satisfy our difference equations we con- 
sider more general correlators than the probabilities of spins. We recover the 
original spin correlators by specializations (see (Q)). We proceed as follows: We 
break the n lines at the chosen row, and change the spectral parameters of the 
lower halves to Then we consider the configuration sum with the 

2n spins fixed to the values ei, . . . , e„ from left to right on the upper row, and 
e'l, . . . ,£'j on the lower (see Figure 1.1). We also consider a 'dislocated' lattice 
obtained as follows. Delete the n lower half lines with the spectral parameters 
Ci I ■ ■ ■ I Cra ■ Therefore, we are left with the n upper half lines whose spectral pa- 
rameters are Ci, . . . , C„. Rename them to ^„+i, . . . , £,2n- Now insert another set 
of n upper half lines with spectral parameters ^1 , . . . , at the next left to those 
which already exist. The orientation of the inserted lines is downward. Thus we 
obtain 2n half lines with the spectral parameters ^1 , . . . , ^2n from left to right 
(see Figure 1.2). We consider the correlator of the 2n spins ti, . . . , r2„ on the 
bottom edges of these half lines, and denote it by F2n{ii, ■ ■ ■ , ^2n)Ti,...,T2„- 
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Figure 1.1 Breaking the successive 
n vertical lines. 



One of our postulates is that the correlator of Figure 1.1 is given by 

In Sect. 1.3 we will give more general statements than this, see (pi]), (p^. If we 
restrict the variables in the above expression as follows 

F^n^\Qn/x, . . . , Cl/X, Cl, . . . , Cn)-e„,...,-e„s„...,.„, (4) 

we obtain the probability of the original n spins taking the values ei, . . . , £„ on 
the regular infinite lattice (see (^). 

Hereafter we will write n, Sj for 2n, ^j, Tj in F^^\ with the understanding 
that if n is odd, the correlator is zero. Set 

i^W(Cl,...,Cn) = E ^^i^HCl,...,Cnk,....s„«s,®---®^^e„ 

£1,...,£ti — ± 

e (g) • • • (8) 

We will denote by Rjk{C) {j < the matrix i?(C) acting on the j-th and the 
k-th. tensor components. We also use the transposition Pjk ~ Rjk{l)- 
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Figure 1.2 The dislocated lattice with 2n 
upper half lines. 

The following are the main results of this paper. In Sect. 1.3 we will give a 
heuristic proof by considering Fn\Ci, ■ ■ ■ Xn)si,...,s„ as a sum over all configu- 
rations with the spins £i, . . . , £„ kept fixed. 

Difference Equation. 

FW(Ci, . . . , xHj, ...,Cn) = Rjj+lix%/Ci+l)-' ■ ■ ■ Rin{x%/Cn)-' 

xRljiCl/Cj) ■ ■ ■ Rj-lj{Cj-l/Cj)Fi'+'\Cl, . . . , 0, • • • , Cn). (5) 

ii-Matrix Symmetry. 

. . , 0+1,0, • • •) = Pij+iRjj+iiQ/Q+i)F^H- ■ ■ , 0. 0+1, • • •)• (6) 

Z2-Invariance and Parity. 

^;l'+'HCl,---,Cnk,...,e„ = ^^«(Cl,...,Cn)-s„...,-e„, (7) 

i^^)(...,-0,...)...,e„... = {-y+'^'e,F^\...,Q,..^...,e,,.... (8) 
Normalization. 

^n+2i- ■ ■ > 0-1, C, ^C, 0' • • ■)---,Sj-i,e,-s,ej,... 

~ Pn\- ■ ■ ' 0-1' 0' • • ■)...,e.i-i,£j,...i (9) 

-^n+2(- ■ ■ ' 0-1' "^C, C' 0' • • ■) ■■■,ej-i,e,£' ,£j 

= 5e,-e,F^\. . .,0-1,0- • • O.-.e^-i.e,,...- (10) 
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1.3 Rotating a half line Let us consider the following more general correla- 
tors; Choose a face of the lattice. Let ni + n2 ~ n (even). Insert ni upper half 
lines and n2 lower half lines into the face, all of which are oriented downward. 
Note that in Sect. 1.2 we considered the cases where ni = n2 or n2 = 0. We 
let the edge on the west side of the said face to be the reference edge used to 
label the sectors, i.e., in the low temperature limit the spin on that edges is 
(— 1)*+-^. (See Figure 2.1.) Denote the correlator of this dislocated lattice by 
F'-'HCi, • ■ • , Cm ; Cl, ■ • ■ , C«2)£i,...,e„i;£i,...,£;2 • Here we place the spectral parame- 
ters and the spins of the upper lines first, from left to right, and then the lower, 
also from left to right. A clue to the difference equations is the following. 

Rotation. 

(Cl 7 ■ • • 7 Crii j Cl , ■ • ■ 7 Cn2)ei,---,£ni ■ ■ .^^,2 
— p( + ) ((j^/x, Cl, . . . , ; ^27 • • • 7 Cn2)-e'i:ei,---,eni \e'2,...,e'^^ , (H) 
= -P"^ (Cl 7 ■ ■ • 7 Cni -1 7 Cl 7 • • ■ 7 Cna 7 Cm / x)ei ,. . . ,e„j^-i;e'j^,. . .,e'^^,-£rii ' (-'-^) 

Let us give heuristic arguments to derive (^)-(|l^). We suppose that the 
correlators are given as configuration sums on the infinite lattice. It makes 
sense to do that because of the normalization of the Boltzmann weights by k. 
We further assume that the vertices on the lattice that are pushed away to 
infinity by a Z-invariant deformation of the lattice, can be neglected. Let us 
use this assumption in a concrete situation. Consider a part of the lattice as 
given in Figure 1.3, on the left. The rest of the lattice, not shown in this figure, 
is irrelevant in the following argument. Denote the configuration sum with the 
spins El and £2 by F{ei,e2', A). The configuration sum for the lattice on the 
right, is given by 

FW(ei,£2;A') = ^ i?(Ci/C2k,e2;ei,4^^^^H£'i7 4; A). (13) 

We can push out the extra vertex in A' to the north as far as we wish by using 
the Yang-Baxter equation, without changing the statistical sum F*^*^ (£1, £2; A'). 
Therefore by the above assumption, we have F*^*) (£1, £2; A') — _F(')(£i, £2; A), 
showing (H). 



7 



Ci Ci Ci C2 





















- 






- 


£1 


/ 

£2 





















- 






- 


£2 


£1 



A 



Figure 1.3 Twisting two vertical lines. 
The extra vertex in A' can be pushed 
up to infinity. 

Next, see Figure 1.4. The right hand side of the figure represents the left 
hand side of (^), the buUet • showing the summation X^e^i- Using unitarity, 
as shown on the left hand side of the figure, we can push away the two vertical 
lines with the spectral parameters ^ and xC,. Therefore, we have (^. 

Eq.(lO) follows from (||) and (||). Eq.(0) is obvious. Eq.(||) follows from the 
parity of the Boltzmann weights. 

Finally, we will show ©-dll). Since the arguments are the same, we will 
consider only (|l|). See Figure 1.5, where C,[ in ( pl| ) is shown as C,. We can 
push away the semi circle with the spectral parameter C,. In the limit, we have 
Eq.dnl). 
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Figure 1.4 The unitarity with the two horizontal hnes 
oriented left, and the case where the orientation of the 
left half of one of the lines are reversed. 





Figure 1.5 The use of the Yang-Baxter equation 
to rotate a half hne. A choice of GSC is written 
exphcitly. Notice that after rotating the GSC changes. 
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2 Spontaneous staggered polarization 

Let us work out the case n = 2 in detail. Since ^2*''(CijC2) depends only 
on the ratio ( = C1/C2, we shall denote F2^'^(Ci,C2) by F^'\c). Set G±(C) = 
F^°\C)±F^^\C)- Eq.(|) for F^"\C) reads G^{x-^C) = ±i?(C)G±(C). Using 

'+(C),G±_(-C) = G|_ 



(0), ® we have Gf_{Q = ±G^+(C), G|_(-C) = G|_(C). Thus the equations 
reduce to 



Set 



G+_(x-^C) = ^(C^)-^ |"f'''l°°^'^;!?'i°° G+_(C). (14) 



f 6 . 4 4\ / 6 . 4 4\ 



{x^z ^;q,x^,x'^)^{qz ^■,q,x'^,x'^)c 



which solves the equation 1^9(2; "^z) — k{z) ^(p{z). It is easy to see that the 
following is a solution of ( p^ ) satisfying the normalization i.e., G^_{x~^) = 
1: 

G+ (C) ^ V?(C^) (a;^a;^)go (g<;g)oo(ga;-iC~^g)oo ^^^^ 

Moreover the solution is unique if we impose further the condition that ^^q < x 
and 

G^_(C) is holomorphic in the neighborhood of x < |C| < x^^ . 

In view of (Q) and following the definition given in Q , the staggared polarization 
Po is 

p^^Gl_{x-^) _{-q-q)Ux^-x^)l, 



G+_(a:-i) (q;g)2^(-x2;x2)^- 

This formula was conjectured by Baxter and Kelland It was proved in the 
trigonometric case g = by Baxter Our result is stronger than this: In the 
context of this section, it says that the ratio of the configuration sums for the 
lattice Aa in Figure 2.1, in the i = sector, with the two different choices of the 
spins at the end edges, i.e., (£1,^2) = (i,T)i is given by 

i g+_(c) + g;_(C) 
2Gt+(C)-G:+(C)' 

and that of the lattice Af, is given by 

i g+_(c/.x) + g;_(c/x) 

2Gt+(C/x)-G:+(C/x)- 
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We have checked these statements to a first few orders in the low temperature 
expansions. 

Aa Aft 

Ci C2 C2 





Figure 2.1 Two types of dislocation. The - shows 
the z = GSC. 



3 Vertex operators 

In this section we shall reformulate the construction of Sect. 1 as an operator the- 
ory. To this end let us first recall the corner transfer matrices (CTMs) following 
. Consider the lattice drawn in Figure 3.1. Number the rows (resp. columns) 
from bottom-to-top (resp. from right-to-left) as —N+ 1, ■ ■ ■ ,0,1, ■ ■ ■ , N. Fixing 
the boundary spins to the i-th GSC, we denote by Ao{(), Ai{(), A2{C), A3{() 
the CTM corresponding to the NE, SE, SW, NW quadrant respectively. When 
formulated in the IRF language they are the transpose of B, A, D, C in p. 
366. For instance the (cr, ct') element of ^3(C) is the partition function of the 
NW quadrant where the horizontal (resp. vertical) boundary spins are fixed to 
a = {ei, ■ ■ ■ ,eN) (resp. a' — {e[, ■ ■ ■ ,e'j^)). From the relations @) it follows 
that ^2(0 = TZAaix-^C^), Ai{C) = 7^A3(C)7^, Ao(C) = A3{x~^C'^)n, where 
TZ — (g) ■ ■ ■ ^ denotes the spin reversal operator. Normalize the CTMs 
so that the largest eigenvalue of A3(C) is 1- Baxter's argument shows that in 
the principal regime we have limjv^oo ^3(C) = C~^j where the operator D is 
independent of C and has discrete eigenvalues 0, 1, 2, • • •. Consequently 

lim Ao iOA, iOA2 (OA, (C) = x^^ . (16) 
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Figure 3.1 The CTMs Aq, 
The 0-th GSC is depicted. 



Denote by Hi the space of eigenvectors of A3(C) in the hmit. In the trigono- 
metric case g = 0, as was shown in |l^, [|], Hi can be identified with the level-1 
highest weight module V{Ki) of the quantized affine algebra U-x{5[2) (with 
~x playing the role of the deformation parameter), and D acts as the grad- 
ing operator in the prinicipal gradation. In the elliptic case, we have no such 
representation theoretical picture. However if Hi ~ ®^^f)HiA is the eigenspace 
decomposition for D, then dmiHi^d cannot change continuously. Hence it should 
be the same as the g = case, giving 



tr^^a;^^ = ^ x^'' dimH»,d = 

d=0 



1 



(a:2;a;4). 



(17) 
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Figure 3.2 Matrix element (<&it^(C))o.cr' 
where a = (ei,- • •,£„), a' = (e'l,- ••,4). 

Consider now the operator whose (cr, u') matrix element is given by Figure 
3.2. We expect that in the large lattice limit ^ oo it will give rise to a 
well-defined operator 

<I>lt'(C) : (18) 

The effect of operating with $*^^(C) amounts to inserting a vertical dislocation 
in the lattice. In view of the expression (^ of the CTMs the function F^^ of 
Sect. 1 is expressed as 
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The arguments in Sect. 1 can be translated to the properties of <i>*^^(Q. We 
summarize them below. 



x^^o$^t'(C)ox-^^ = $^r(.^C), (20) 

'^l+leAC2)^rtliCl) = E ^(Cl/C2)e.s.;e;si<J>:+le;(Cl)1>:^(C2), (21) 

There exists a linear isomorphism u : TYq — ^ 'Hi such that 

,yoD = Do,y, $1,(0 =z^o$0_^(C)oj., (22) 

C(-c) = (-!)'+'£ $:r(c), (23) 

E '^r+ieicWit\c) = id, $^t^(c)$:;ie'(c) = s,,, x id. (24) 

e=± 

In eq.(|2|) we have set $*^"'"i(C) = ^'f^x~^C,). The difference equations (5h for 
( p^ is an immediate consequence of (pO[), (^l|). In the trigonometric case (1^) is 
related to the vertex operators in |l|, | via ^l+\0 = g-'^/^C^'^+^y^^'^^^^^C) 
where g is a scalar independent of C. 
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